UMD-PP-03-047 



Lepton Flavor Violation and Neutrino Mixings in a 

3x2 Seesaw Model 

Bhaskar Dutta 1 and R.N. Mohapatra 2 

1 Department of Physics, University of Regina, Regina, SK, S4S 0A2, Canada 
2 Department of Physics, University of Maryland, College Park, MD 20742, USA 

(Dated: May, 2003) 

Abstract 

We study the predictions for BR(u — ► e+7) and BR(t — > e, /i+7) in a class of recently proposed 
horizontal SU (2) model that leads to a 3x 2 seesaw model for neutrino masses. We consider two such 
models and obtain the correct bi-large mixing pattern for neutrinos. In these models, the effective 
low energy theory below the SU(2)h scale is the MSSM. Assuming a supersymmetry breaking 
pattern as in the minimal SUGRA models, we find that consistent with present g-2, b — ► s + 7 and 
WMAP dark matter constraints on the parameters of the model, the BR{u — » e + 7) prediction 
is accessible to the proposed MEG experiment at PSI making this class of models testable in the 
near future. 
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I. INTRODUCTION 



Evidences for neutrino oscillations and hence for neutrino masses and mixings seem to be 
getting more and more solid. The results from solar and atmospheric neutrino experiments 
as well as those using terrestrial neutrinos as in the K2K and KAMLAND experiments 
which are responsible for this very important conclusion have together with CHOOZ and 
PALO- VERDE experiments clearly established the mixing pattern among different neutrino 
generations. It seems that two of the three neutrino mixings i.e. 9 12 and #23 are large and 
the third angle #13 is smallQ]. We will follow the literature and call this the bi-large mixing 
pattern. 

This good news on the mixing front is however not shared by the masses i.e. our knowledge 
of the neutrino mass pattern is far from clear. At the moment three obvious possibilities i.e. 
(i) normal mass hierarchy, (ii) inverted mass hierarchy and (iii) degenerate mass pattern are 
equally viable. The challenge for future experiments is therefore twofold: to pin down the 
mass pattern more definitively and to improve measurements of the three angles to more 
precise values so that different theoretical models can be tested. For example it has been 
noted that if the present upper limit on the angle #13 drops by one order of magnitude, a 
lot of very interesting models for neutrino masses will be ruled out. Of course more models 
being ruled out makes the direction of physics beyond the standard model that much more 
clear. It is therefore important to extract the predictions of the various models as far as 
possible to narrow the choices. 

One of the very intriguing mass patterns for neutrinos is the inverted pattern, where the 
two heavier neutrinos are split by the solar mass difference Am^, whereas their absolute 
mass value is the square root of the atmospheric mass difference square i.e. Am\. In the 
extreme limit where Am?, = 0, this mass matrix admits an L e — — L T symmetry. If 
this pattern is confirmed by experiments, the existence of the L e — L M — L T symmetry could 
provide a clue to new local symmetries that operate at the higher scales where neutrino 
masses arise. 

In two recent papers)^, it was shown that if there is a leptonic SU(2)h horizontal 
symmetry^] acting both on the left as well as the right handed leptons of the standard 
model, then the theory must have two right handed neutrinos to cancel the global SU(2)h 
Witten anomaly. This provides a rationale for the existence of the right handed neutrino 
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needed for implementing the seesaw mechanism. This leads to a 3 x 2 seesaw formula for 
neutrino masses. The resulting neutrino mass matrix naturally leads to an inverted mass 

n 

pattern for the light neutrinos after seesaw mechanism 4 and in a certain areas of the pa- 

n 

rameter space gives rise to an approximate L e — — L T symmetry |5|] of the neutrino mass 
matrix. Thus if the inverted mass pattern for neutrinos is confirmed by experiments, the 
SU(2)h models must be given very serious consideration. 

As discussed in ref . 3| , there are two possible patterns for the charged lepton mass matrix 
in these models and since the PMNS mixing matrix gets contributions both from the charged 
lepton as well as the neutrino sector, there will be two sets of predictions for the neutrino 
mixings. In this paper, we focus on them separately and study their physical implications, 
especially their predictions for lepton flavor violating branching ratios BR(fi — > e + 7) and 
BR(t -> e,/x + 7). 

The key point is that below the horizontal symmetry scale, the theory is the minimal 
super symmetric standard model (MSSM) and as is well known Jf], Q, 0, 0], the presence of 
lepton mixings leads to mixings between sleptons via the renormalization group evolution 
even though at the GUT scale all slepton masses being equal forbid any such mixings. These 
in turn lead to lepton flavor violating decays such as BR(fi — > e + 7) and BR(t — > e, \i + 7) 
through one loop graphs involving the gauginos. 

In order to arrive at our result, we first fix the model parameters such that they lead 
to observed neutrino mixings. We then find that for both the models referred to above, 
in most of the allowed parameter space for MSSM, the BR(\i — > e + 7) is around 10~ 12 , 
which is within the range of the next generation of experiments |lfl] . The constraints on the 
supersymmetry parameter space come not only from the collider searches for supersymmetry 
and the LEP bound on Higgs boson mass but also from the Brookhaven g-2 results as well 
as the WMAP dark matter constraints [llj|. The last constraint arises from the lightest 
neutralino Xi as the dark matter candidate and using the recent WMAP results for the dark 
matter content of the universe i.e. 0.094 < VLh 2 < 0.129 (2a). In deriving these results, 
some assumptions are needed for supersymmetry breaking parameters of the model. We 



use the mSUGRA choice[12j at the GUT scale to obtain the constraints i.e. universal mass 
for superpartners, common gaugino mass and A-terms being proportional to the Yukawa 
couplings. It is important that the GUT scale is higher than the seesaw scale to get sizable 
flavor mixing effects. 
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This paper is organized as follows: in sec. 2, we briefly review the salient features of the 
SU(2)h models and write down the 3x2 seesaw formulae for the neutrino masses and mixings 
as well as the charged lepton mass matrices. We also present the exact neutrino mixing 
matrices in the limit of zero electron masses. In sec. 3, we present the numerical analysis 
of the neutrino mixings and masses and write down the choices of the model parameters 
that are phenomenologically acceptable. In sec. 4, we discuss the predictions for lepton 
flavor violating branching ratios in the model for these choices of parameters. In sec. 5, we 
conclude with a summary of the results. 

II. MODELS FOR LEPTON MIXING 

We consider extensions of the MSSM where the gauge group extended to SU(3) C x 
SU(2)l x U{1)y x SU(2)h above the seesaw scale. Below the seesaw scale the model is 
MSSM so that the usual results for gauge coupling unification and electrweak symmetry 
breaking hold. 

The basic motivation for considering these models is to find a symmetry rationale for 
the existence of right handed neutrino othet than the fact that it is needed to implement 
the seesaw mechanism. An elegant extension of the standard model that brings in the 
right handed neutrino automatically is to include the local B-L symmetry in the theory i.e. 
theories based on the gauge group £77(2) l x U (1)i 3R x U(1)b-l ot SU (2) L x SU (2) R xU(l) b-l 
or the GUT group SO(10). Since our understanding of neutrino physics is not complete, it 
is important to explore other possibilities. 

A different possibility is to look for horizontal (or family) symmetries that arise in the 
standard model in the limit of vanishing fermion masses. To see the possibilities, let us work 
with the extension of the standard model based on the gauge group Gstd x Ghi where Gh 
actos on the different generations. Three interesting possibilities that provide independent 
motivations for the existence of right handed neutrinos are: (i) Gh = S77(3) 9+ ^|]J3]; (ii) 
Gh = SU(2) q+ £ |2jand (iii) Gh = SU(2) e ^ with the obvious choice of fermion assignments. 
Here we will work with the case (iii) where we assume that the quark superfields transform 
under the SU(2) H group as siuglets aud leptons a s doublets. In this case, it was noted in 
|2| that while this model is free of local anomalies, absence of global Witten anomalies [14] 
implies that one must add two right handed neutrinos to the theory. Thus the existence 
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of the right handed neutrinos are required for the consistency of the theory and not just 
to implement the seesaw mechanism. This is the model we want to study in detail in this 
paper. 

The assignment of the leptons and Higgs superfields under the gauge group SU(3) C x 
SU{2) L x U(1) Y x SU(2) H is given as follows: 



Table ] 




= (L e , Lfj) 


(1,2,-1,2) 




(1,2,-1,1) 


E c = (fi c , -e c ) 


(1,1,-2, 2) 




(1,1,-2, 1) 


N c = (^,-u c e ) 


(1,1,0,2) 


K 


(1,1,0,1) 


Xh = ( Xi Xi ) 


(1, 1, 0, 2) 


Xh = ( -X2 Xi ) 


(1,1,0,2) 


H u 


(1,2,1,1) 


H d 


(1,2,-1,1) 


A H 


(1,1,0,3) 



Table caption: We display the quantum number of the matter and Higgs superfields of 
our model. 

Here L e ^ T denote the left handed lepton doublet superfields. We arrange the Higgs 
potential in such a way that the SU(2)h symmetry is broken by < xi >=< X2 >= vhi] < 
X2 >=< Xi >= vh2 and < >= v' H , where vn,i,v' H ^> v w k- The vevs for x are chosen 
so as to cancel the D-terms and leave supersymmetry unbroken below the scale vh- Note 
that we have used the SU(2) H symmetry to align the A H vev along the Ih,3 direction. At 
the weak scale, all the neutral components of the fields H u and Hd acquire nonzero vev's 
and break the standard model symmetry down to SU(3) C x U(l) em . We denote these vev's 
as follows: < >= kq and < >= K cot(3; Clearly k has values in the 100 GeV range 
depending on the value of tan(3. 

Note that < A H >^ breaks the SU{2) H group down to the C/(l)i e _r group which is 
further broken down by the xh vevs. Since the renormalizable Yukawa interactions do not 
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involve the xh field, this symmetry (L e — L M ) is also reflected in the right handed neutrino 
mass matrix. This turns out to play a crucial role in leading to a bimaximal mixing pattern. 

To discuss the detailed leptonic mixings, we now write down the Yukawa superpotential 
of the theory. The two models we discuss arise from the two ways the horizontal doublets 
couple to fermion superfields. In one case, we will allow only the \h field to couple to all 
fermions (to be called model I) and in the second case, we will allow the xh to couple to 
the right handed neutrino sector only and xh to couple to the charged leptons only. This 
will be called model II. 



A. Model I 



The Yukawa superpotential for this model is given by: 
W Y = ho^H^ + L^H^+hL^u^ + KxiWM-zfN^r-AnN' 



(1) 



^(L eX 2 - L, X i)H d T c + ^L T H d (fi c X 2 + e c X i) + h' 3 L T H d r c + h' 2 (L e e c + L^ c )H d 

This can be made natural by a choice of the discrete Z 2 symmetry under which x~h is odd 
and all other fields even, so that it does not couple to matter fields in the superpotential. < 
A° H >= v' H directly leads to the L e — L M invariant v^—v^r mass matrix at the seesaw scale, as 
already noted. The xh vev contributes to this mass matrix only through nonrenormalizable 
operators and we assume those contributions to be negligible. We also do not include 
any term where A H couples to light fields. Since the theory is supersymmetric, any term 
omitted from the superpotential will not be induced by loops due to the non-renormalization 
theorem. Similarly there will also be some small contributions from the v t r sector if we did 
not decouple it completely. We ignore these contributions in our analysis. Further, we define 
« lj2 = <Xl * >Ko -1-100 GeV. 

To study neutrino mixings in this model, let us write down the 5x5 mass matrix for 
neutrinos: 



M, 












h K 

















h K 











h\Ki 


hiK 2 


h K 





hiKi 





fV H 





h K 


hiK 2 


fv' H 






\ 



(2) 



HQ 



This is a 3 x 2 seesaw matrix! 2, 15] and after seesaw diagonalization, it leads to a 3 x 3 light 
neutrino mass matrix as follows: 



M v = -MoM^Ml 



(3) 



where Mp 



h K 
h n Q 



i 



y hiKt h l K 2 j 

mass matrix M u is given by: 



1 

1 



. The resulting light Majorana neutrino 



M v 



1 

7^ 



(4) 



(h K ) 2 h hiK K 2 ^ 
(/i /t ) 2 hoh^o^ 
i h hiKoK,2 h hiK K,i 2h\K\K 2 j 

To get the physical neutrino mixings, we also need the charged lepton mass matrix defined 
by iphMnipR. This is where the difference between the two models appear. 
In model I, the charged lepton mass matrix is given by 

' h' 2 Ko —h! x K2 ^ 

h' 2 Ko h' x Ki 



Mi 



cot {3 



(5) 



y h'^Ki h' 4 K 2 h' 3 K J 

In order to study physical neutrino mixings, we must diagonalize the M v and Me matrices. 
In ref.|2|, we diagonalized the neutrino mass matrix analytically in the limit m e = 0. Let us 
discuss it for completeness. Since h 2 K Q is proportional to the electron mass, we will ignore 
temporarily. Then, defining the matrices that diagonalize the charged lepton mass matrix 
as D e = U ( e L) M e U ( e R) \ we get 

( si ci o N 

CpC\ -CpS! Sp 



u, 



(6) 



where two angles 612 are given by: 



sin9\ = sj. 



Ki 



(7) 



sin9 2 = s 2 



h K 
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m M /m r and we have ignored terms of order m e jm ll . Similarly the 



and tan2/3 ~ 

matrix that diagonalizes M. v is given by 

^ C\C' — SiS 2 s' C\S' + SiS 2 c' C 2 Si ^ 

—S\d — C1S2S' —sis' + C1S2C 1 c 2 c\ 



V 



(8) 



-c 2 s' c 2 d -s 2 / 



where s' = sin9' with 6 1 ' given by tan26' = s ^ ■ The PMNS matrix that is measured 

in neutrino oscillation is given by Tj\ U v and has the form, in the limit of m e = 0: 



U 



PMNS 



V 



(9) 



s 2 s s 2 c c 2 

-(dCp + C 2 s'sp) -{Ops' - SpC 2 d) -S 2 Sp 

'ds p - e 2 s'cp) (c 2 dcf3 + s'sp) s 2 cp J 
Using this, it was concluded in 2J that there is an approximate relation between the neutrino 
parameters given by: 



U^cos26 Q 



2Am 2 A 



(10) 



It is clear from the above equation that for a given value of Am^/Am^, the smallest value 
of U e s will happen for the largest value of sin 2 26 & . In this paper first we numerically solve 
for the neutrino mixings and then proceed to calculate the lepton flavor violation effects for 
the allowed choices of the parameters of the model. 



B. Model II 



In this case, we impose a discrete symmetry under which the fields (e c , /i c , r c ) and xh are 
odd and all othee fields are even. The Yukawa superpotential invariant under this is: 

Wy = hoiL.H^ + L^H^ + hMu^ + ^X^HjM-tfN^T.r-AnN^ (11) 

jj(L e Xi + L,X2)H d T c + ^L T H d (fi c xi ~ e c x 2 ) + h' 3 L T H d r c + h' 2 (L e e c + L^ c )H d 

Note that the vev's < Xi,2 >= ^2,1 for the D-terms to cancel at the high scale. Substituting 
these vevs and the vev's of H u ^ d , we get exactly the same neutrino mass matrix as in model 
I i.e. Eq0]but the charged lepton mass matrix becomes: 

^ h' 2 KQ h[Ki ^ 
M e = cot(3 h' 2 K h[K 2 ■ (12) 
y -h' 4 n 2 h' 4 Ki h' 3 K j 
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The matrices that diagonalize this matrix are different from the first case (model I) and the 
Jjf^ in this case is given by: 



U, 



(L) 



Ci —Si 

\ SpSx SpCi C/3 J 



The Upmns in this case is given by 

/ 

UpMNS — 



d s' 

s's 2 Cp - C 2 s'sp S 2 c'Cf3 + C 2 c'sp C 2 C/3 ~ S 2 Sp 

\ {-s's 2 sp - c 2 s'cp) {c 2 dcp - s 2 c'sp) S 2 Cp - C 2 Sp ) 



(13) 



(14) 



Note that we have neglected corrections of order ^m e /m^ and once they are turned on, we 
indeed have a nonzero C/ e 3 and small. We will discuss the detailed numerical predictions in 
the next section. 



III. NEUTRINO MASS FIT 



A. MODEL I 



The expressions for the Dirac neutrino matrix and the charged lepton matrix are given 
in the previous section. We will use those expressions in order to derive the light neutrino 
masses and the mixing angles numerically. That will fix the model parameters that will go 
into the renormalization group evolution of the slepton masses. We will generate the fit for 
different values of tan/3 and model parameters. For example, at tan/3 = 30, the charged 
lepton mass matrix is given by (at the scale where the horizontal symmetry gets broken) : 

/7.5xl(T 3 -1.135N 



(15) 



M t = 7.5 x 1(T 3 0.247 

v 2.14xl0- 2 0.0984 0.85 / 

This matrix needs h[ = 1.135 xl0 _1 , h' 2 = 1.29 x 10~ 3 , h' 3 = 1.466 x 10" 1 and h' A = 9.8 x 10~ 3 
with k 2 = 10 and ^ = 4.588. We get the correct values of charged lepton masses at the 
weak scale from the above matrix. We use the MSSM RGEs between the horizontal scale 
and the weak scale. Using the same ks, ho = 9.2 x 10~ 2 and hi = —2.92 x 10" 1 , we find the 
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light neutrino mass matrix (in eV): 



My 



5.15 x 1(T 2 -9.3 x 1(T 3 \ 

5.15 x 10~ 2 -2.04 x 10~ 3 . 

\-9.3xl0- 3 2.04 xl0~ 3 7.4 x 10~ 4 / 

The PMNS matrix and neutrino eigenvalues at the weak scale are : 

/ 0.806 0.58 -0.109 \ 
Upmns = -0.27 0.527 0.805 . 

V 0.526 -0.619 0.582 J 

(mi, m 2 , m 3 ) = (5.276 x 10~ 2 , 5.202 x 10~ 2 , 0) eV . 



(16) 



(17) 



(18) 



From the above fit, we find sin 2 26 Q = 0.87, Am% = 7.7 x 10~ 5 eV 2 and sin 2 26 A = 0.88, 
Am\ = 2.7 x 10~ 3 eV 2 which are well allowed by the experimental results. 

The above fit is obtained for a particular value of We can change that value and the 
Yukawa couplings hs and h's and still can obtain very good fits. In fig.l we show the model 
points which are within the experimental limits of sin 2 26 & and sin 2 26a- We find many model 
points which satisfy the allowed ranges. In Fig.2, we show the \U e3 \ as a function of sin 2 26 Q 
for the model points which satisfy the experimental limit on sin 2 26a- The range of U e3 is 
0.07-0.13 in this model and is allowed by the current experimental upper bound 0.22|l6j. In 
fact, it is expected that the planned long base line experiments including MINOS have the 
ability to probe this range in the next few years. 

The scale of the horizontal symmetry breaking is assumed to be 3 x 10 13 and the magnitude 
of the Majorana coupling is / ~ 1.6 x 10 _1 . The fit does not change if we move the scale 
up or down. 

The fit for tan/3 =10 and 40 look similar. The charged lepton mass matrices for these 
two tan (3 are given below. 



For tan f3 = 40 we have: 



Me 



( 8.2 x 10~ 3 


V2.33 x 10~ 2 





8.2 x 10~ 3 
1.07 x 10- 1 



-1.30 
2.83 x lO" 1 
9.5 x 10- 1 



(19) 



We keep the light neutrino matrices unchanged and we find sin 2 26 Q = 0.88, Am^ = 7.7 x 
10~ 5 eV 2 and sin 2 26 A = 0.87, Am 2 A = 2.7 x lO" 3 eV 2 . 
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1 

0.95 
0.9 
0.85 



8 0.8 

d 



05 



0.75 
0.7 

0.65 
0.6 



0.6 0.7 



0.8 

sin 2 20 




k** *** 



**** ** 



0.9 



atmos 



FIG. 1: sin 2 29Q vs sin 2 28 atmos - The stars show the model points. The horizontal and vertical 
solid lines show the current experimental limits. 



For tan (3 = 10 we have: 



Mi 



( 6.8 x 1(T 3 

6.8 x 1(T 3 

\2.088 xl(r 2 9.58 x 1(T 2 



-9.9 x lO- 1 ^ 
2.15 x KT 1 
8.0 x 10- 1 / 



We keep the light neutrino matrices unchanged and we find sin 2 29 & = 0.88, AttIq 



(20) 



7.7 x 



10" 5 eV 2 and sin 2 26, 



0.92, Am 2 A 



2.7 x 10" 3 eV 2 



B. MODEL II 

In this model, the Dirac neutrino matrix is similar to what we have in the previous model. 
The charged lepton matrix is the transpose of that in the other case. We can again obtain 
the neutrino fit for different values of tan (5. For example, at tan (3 = 30, The charged lepton 
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0.95 
0.9 
0.85 



8 0.8 

d 



05 



0.75 
0.7 

0.65 
0.6 



******* 

* # 

***** * 
**-«***■* 
******* 
*>***** 

***-***** 



0.1 0.2 0.3 0.4 

IU J 

e3 



FIG. 2: sin 2 29 & vs | C7" e 3 1 . The stars show the model points. The horizontal solid lines show the 
current experimental limits. 



mass matrix is given by (at the scale where the horizontal symmetry breaks down): 

/8.28xl0" 2 5.0 xl(r 5 \ 

M e = 8.28 x 10~ 2 4.0 x KT 3 . (21) 

V -1.44 0.018 6.9 x 10- 3 j 

This matrix requires ^=0.0143, h' 3 = 1.2 xl0~ 3 , h[ = l.OOxlO" 4 , h' 4 = 3.6xl0~ 2 for k 2 = 40 
and ^ = 0.5. The charged lepton mass matrix produces the correct values of charged lepton 
masses at the weak scale after using the RGEs. Using the same k's, ho = 7.67 x 10~ 2 and 
hi = 3.07 x 10 _1 we find the light neutrino mass matrix: 



3.27 x 10~ 2 \ 
4.09 x 10~ 4 



/ 3.55 x 10~ 2 

3.55 x 10" 2 
\3.27 x lO" 2 4.09 x KT 4 7.5 x 10" 4 J 
The PMNS matrix and neutrino eigenvalues are : 

/ -0.724 0.68 



(22) 



U 



PMNS 



0.527 0.513 



-0.042 \ 
-0.677 



(23) 



y -0.44 -0.512 -0.734/ 
12 



(mi, m 2 , m 3 ) = (4.87 x 1(T 2 , 4.79 x 10~ 2 , 0) eV . (24) 

The solar angle prediction i.e. sin 2 26 Q = 0.95 for this model is on the higher side and there- 
fore as the experimental error on 612 narrows, this model will be tested. The atmospheric 
mixing angle sit^IQa is 0.98 and is well allowed by the experimental results. The prediction 
for U e3 = —0.042 and it will also be probed in the near future. 

In the following, the scale of horizontal symmetry breaking our calculations is assumed 
to be 3 x 10 13 and the magnitude of the Majorana coupling is / ~ 1.6 x 10 _1 . 



IV. LEPTON FLAVOR VIOLATION 



It is well known that in the standard model where neutrino masses vanish, there is no 
lepton flavor violation. Once the neutrinos have masses as well as mixings, there are small 
amounts of lepton flavor violation leading to processes such as £j — > + 7 with i 7^ j. 
However, if only a Dirac mass( by adding a right handed neutrino) or a Majorana mass (by 
adding a Higgs triplet) of the neutrino are added to the standard model, then lepton flavor 
violating processes which always arise at one loop or higher loop level in gauge theories 
have branching ratios proportional to ( ^f- ) and are therefore very tiny. In models with 

\ m w ) 

supersymmetry however, the situation can change drastically. In the presence of neutrino 
masses, the superpartners of leptons could in general mix and the LFV branching ratio in 

/ ml \ 2 

that case becomes of order \—Jh J where m^r is a typical mass of the slepton and since the 
slepton masses are expected to be of order of the weak scale, the LFV branching ratios in 
general can be quite large. The actual magnitude however depends on the specific model 
for neutrino masses as well as on the pattern of supersymmetry breaking. If we assume the 
neutrino masses to arise from the seesaw mechanism, then depending on how the seesaw 
mechanism is implemented in a theory, the predictions for LFV can be different. 

Similarly, the nature of supersymmetry breaking is also important since in the super- 
symmetric limit the radiative flavor changing decays vanish. To illustrate this point, let us 
consider a simple seesaw model where the susy breaking slepton masses are universal at the 
seesaw scale and the SUSY breaking A terms are proportional to Yukawa couplings as in the 
minimal SUGRA models. We can go to a basis where the charged lepton Yukawa couplings 
are diagonal. This simply rotates the neutrino mass matrix and keeps all other terms in the 
Lagrangian flavor In such a theory, just below the seesaw scale, flavor mixing is present only 
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in the neutrino masses from the seesaw formula. As we extrapolate the slepton masses down 
to the seesaw scale, they remain diagonal since the charged lepton Yukawa coupling is diag- 
onal. The only source of lepton flavor violation is then in the neutrino masses. Therefore, in 
this case BR(fi — > e + 7) oc ( -^f- ) and is therefore negligible. To obtain a significant LFV 

\ m w J 

effect therefore, one must assume that the universality of slepton masses and the propor- 
tionality of the A terms must hold much above the seesaw scale. Fortunately, this is what 
is generally assumed in minimal SUGRA models anyway. Thus no observable LFV effect 
would simply be a manifestation of the fact that it is the seesaw scale where universality of 
slepton masses holds e.g. in models like gauge mediated supersymmetry breaking. In this 
paper, we will assume that the universality of scalar masses holds at a scale higher than the 
seesaw scale. 

Turning to our case, since we only have an SU(2) horizontal symmetry, we can dispense 
with any assumption regarding universality and allow for the most general SUSY breaking 
pattern consistent with the SU(2)h symmetry. For definiteness, we will assume universality 
of scalar masses at a higher scale (e.g. the GUT scale) and use RGE effects of the lepton 
Yukawa couplings to get the slapton masses and the A terms at the vh scale. This is a non- 
universal profile at the seesaw scale and leads to function of the Yukawa couplings. 
If we then rotate i.e. m~ = m~ = m 2 ~ m 2 , + ( 6+2a ^o+ fe 2 ) m o £ n A% anc i m 2 _ m 2 _|_ 
(g+ggjjftg )™ g n My__ we nQW ro ^ a t e the charged lepton Yukawas just below the horizontal 
symmetry scale to go to a diagonal basis for charged lepton Yukawa couplings, we will 
generate slepton flavor mixings such as rn 2 LLlz and Tn 2 LL23 - A combination of these two 
terms at the weak scale will give rise to a nonvanishing branching ratio for fi — > e + 7 and 
r — > e, n + 7. If we focus the slepton masses alone, we roughly get 



B0i->e + 7)« ™ L_L* &l (25) 



One can also in a similar manner estimate the contribution of the A terms. In fact it turns 
out that it is the A-term slepton mixings that make the dominant contribution to the LFV 
branching ratios in our model. 

To proceed with the calculations for LFV processes in our case, we will use mSUGRA 
models which has five parameters at the GUT scale i.e. mo, mi/2, Aq, tan/3 and yu[12]. We 
will assume that at the GUT scale all scalar masses have the common value m as noted 
in the previous section; gauginos have the common mass mi/%] the trilinear SUSY breaking 
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terms are proportional to the Yukawa couplings and the rest of the theory is defined by the 
the superpotentials in section II. First note that the renormalization group evolution (RGE) 
is caused only by dimension four terms in the Lagrangian. The charged lepton Yukawa 
matrix responsible for the RGE from the GUT scale to the weak scale is diagonal since the 
only dimension four Yukawa couplings are h' 2 3 . Similarly, the dimension 4 Dirac neutrino 
couplings are also flavor diagonal with the Yukawa coupling ho. Thus as we run down from 
the GUT scale to the horizontal symmetry breaking scale vh, the slepton masses remain 
horizontal; The A terms get modified. 

The RGE from the vh scale down to the weak scale is the MSSM RGE. We generate 
the off diagonal elements of the charged lepton Yukawa matrix at the horizontal scale from 
the nonrenormalizable operators and Dirac neutrino matrix becomes a 3 x 2 matrix. The 
off diagonal elements of the charged lepton mass matrix induce the flavor violating terms in 
the m 2 LL , and m RR mass matrices via the RGEs from the horizontal scale down to the 
weak scale. For example, E Al2 term (which contributes to m 2 LR ) gets generated via terms 
E A U}U h UiU}E A pieces in the RGE (^oc ^{AE A U}Ui + WiU}E A ), where U t is the 3 x 3 
matrix for the Dirac Yukawa coupling). We have £7i3,3i and t/32,23 terms in both models. 
These terms generate large BR(/i — > e + 7) . 

The calculation of BR(jU — > e + 7)involves both neutralino and chargino diagram contri- 
bution. In some parameter space one type dominates and the other type dominates in some 
other. Before we discuss the BRs. of the lepton flavor violating decay modes, we need to 
discuss the supersymmetric parameter space. 

A. Supersymmetry parameter space 

The mSUGRA parameter space is constrained by the experimental lower limit on nih and 
measurements of b —>■ 57 and recent results on dark matter relic density [h,]. For low tan/3, 
the parameter space has lower bound on mi/2 from the light Higgs mass bound of rrih > 114 
GeV. For larger tan (5 the lower bound on mu2 is produced by the CLEO constraint on the 
BR(6 — > S7). The recent WMAP results have led to stronger constraints on dark matter 
density which in turn reduces the allowed parameter space of MSSM mostly to the co- 
annihilation region for mo, < 1000 GeV. As is well known, co-annihilation requires 
some other superpartner to have mass close to the lightest neutralino. In the mSUGRA 
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model co-annihilation happens between the lightest slepton and the lightest neutralino 17]. 
One can see this analytically for low and intermediate tan f3 for mSUGRA where the RGE 
can be solved analytically ,18]. At the electroweak scale one has for the average slepton 
masses 

rhln = ™o + (6/5)/im? /2 - sin 2 6 W M^ cos(2/3) (26) 

m^o = (a 1 /a G )m 1 / 2 (27) 

where /* = [1 - (1 + A*)~ 2 ]/A, t = \n(M G /M z ) 2 and (3 X is the £7(1) (3 function. Numerically 
this gives for e.g. tan (3 = 5 

m\ R = m 2 + 0.15m 2 /2 + (37GeV) 2 

m|o = 0.16m 2 /2 (28) 

Thus for m = 0, the e# becomes degenerate with the Xi a t m i/2 = 370 GeV, i.e. coanni- 
hilation effects begin at m 1( / 2 = (350 — 400) GeV. For larger m^, the near degeneracy is 
maintained by increasing mo, and there is a corridor in the mo — mi/2 plane allowing for an 
adequate relic density. For larger tan/3, the stau is the lightest slepton and its mass is close 
to the neutralino mass for a large region of parameter space. 

The previous relic density bound 0.07 < Qh 2 < 0.21 reduced the parameter space mostly 
into coannihilation bands. The most recent WMAP result 0.094 < Qh 2 < 0.129 has narrowed 
the co-annihilation regions further Iil9ll . 

n 

Further constraints on the parameter space arise from the present muon g-2 results [2PJ, if 
we take the deviation from the standard model value as an upper limit on the supersymmetric 
contribution to g-2. For instance, for tan (3=10, we get m 1( / 2 < 430 GeV. Since the Higgs mass 
implements a lower bound rrii/2 >300 GeV, we do not have much parameter space left for 
this value of tan (3. However, for tan (3=30 and 40, the upper bounds are < 720 and 790 
GeV respectively from the g-2 of muon data and thus the allowed parameter space becomes 
bigger. We will therefore restrict ourselves to larger values of tan (3 for the discussions of 
flavor violation BRs. 

In Figs. 3 and 4, we show the BR(/i — > e + 7) for tan /3=30 and 40 for different values of 
Aq in model I. We tune mo such so that Qh 2 constraint is satisfied. This mostly happens 
where the light stau mass is very close to the lightest neutralino mass. For example, for 
m i/2 — 400 GeV, m ~ 208-215 GeV satisfies the Qh 2 constraint for tan/3 = 40 and A = 0. 
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FIG. 3: BR[/i — > e + 7] as a function of m.1/2 f° r different values of ylo f° r model I. 

The BR remains almost unchanged over the small range of mo. The parameter space shown 
in the figures are allowed by the BR[6 — > 57] constraint. From these figures, we find that as 
mi/2 increases BR decreases. In the figures, the dependence of the BR on A$ is also shown. 
The relative sizes of the chargino and neutralino diagram contributions change with A . The 
BR, in most of the parameter space, is around 10~ 12 . Choosing a smaller ^a. ratio dcreases 
the BR ratio further. But this ratio has a lower bound from the fit to the solar mixing angle. 
The lowest possible BR is a factor of 2 smaller than what is shown in the figure. The BR 
also can be reduced by a factor of 4 if we increase the seesaw scale from 3 x 10 13 to 10 15 

Turning to the BR(r — > //(e) +7), it is at least 3 order magnitude below the current 
experimental result in the parameter space allowed by the BR(/i — > e + 7). 

In Figs 5 and 6, we show the BRs. for fi — > e + 7 for tan (5=30 and 40 for A$ =0 and 
A = 200 GeV in model II. The mo is chosen such so that flh 2 constraint is satisfied. The 
t — » /x(e) + 7 BR is at least 3-4 order magnitude below the current experimental limits. 

Furthermore, the parameter space where the BR(/i — > e^y) is larger, the other observables 
e.g. BK(B S — > fifi) can be observed in the current RUN at the Tevatron. Similarly, for these 
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FIG. 4: BR[/u — > e + 7] as afunction of TO1/2 for different values of Aq for model I. 



values, cr^o_ p are also large |21| and can therefore be observable in planned dark matter 
experiments. Thus, once the parameter space is narrowed from these different measurements, 
this model can be deciphered by the /1 — > e + 7 experiment. 



V. CONCLUSION 



In summary, we have analysed two models with local SU(2) horizontal symmetry that lead 
naturally to bi-large mixing pattern for neutrinos. We determined the parameters of these 
models from the detailed numerical fit to the neutrino masses and the mixing angles. Using 
these models we then obtained predictions for lepton flavor violating decays BR[/z — > e + 7] 
in mSUGRA framework. The parmeter space, we consider, is constrained by the following 
experimental constraints: BR[6 — > sj], Higgs mass and the recent relic density results from 
WMAP and the Brookheven g-2 experiment. We find that, in both models, the branching 
ratio is large in a large region of parameter space and is in the range accessible to the current 
round of searches for this process. This parameter space can simultaneously be explored at 
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FIG. 5: BR[/i — > e + 7] as afunction of roiw for different values of j4q for model II. 



Tevatron, LHC and the upcoming of dark matter detectors. 
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